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Abstract 


Two  classes  of  finite  and  infinite  moving  average  sequences  of  bivariate 
random  vectors  are  considered.  The  first  class  has  bivariate  exponential  mar¬ 
ginals  while  the  second  class  has  bivariate  geometric  marginals.  The  theory  of 
positive  dependence  is  used  to  show  that  in  various  cases  the  two  classes  consist 
of  associated  random  variables.  Association  is  then  applied  to  establish  moment 


inequalities  and  to  obtain  approximations  to  some  joint  probabilities  of  the 


related  bivariate  point  processes, 
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1.  Introduction  and  Summary. 

A  reigning  stationary  model  in  time  series  analysis  is  the  p*l  moving  average 
(M.A.)  model  given  by: 

(1.1)  X(n)  -  £°°  A(j )  e  (n-j),  n-  0,  +1, +2,..., 

j— 00 

where  A(j),  j *  0,+l,+2,. . . ,  is  a  sequence  of  p*p  parameter  matrices  s.t. 

I  II  A(j )||  <  00 ,  and  e(n),  n- 0,+l,+2,. . .  is  a  sequence  of  uncorrelated  p*l 

j— 00 

random  vectors  (r.ve.s)  with  mean  zero  and  common  covariance  matrix.  It  is  well 
known  that  this  model  emerges  from  many  physically  relizable  systems  (see  for 
example  Hannan  (1970),  p.  9),  however  in  some  physical  situations  where  the 
r.ve.s  X(n)  are  either  positive  or  discrete  the  preceeding  assumptions  about  the 
r.ve.s  e(n)  are  Inappropriate  (see  Lewis  (1980),  p.  152). 

Several  researchers,  addressing  themselves  to  this  problem,  have  been  con¬ 
structing  univariate  stationary  M.A.  models  and  univariate  stationary  auto¬ 
regressive  moving  average  (A.R.M.A.)  models  where  the  r.v.s  X(n)  have  exponential 
(exp. )  or  gamma  distributions  (dist^s)  and  discrete  models  where  X(n)  assumes 
values  in  a  common  set.  Lawrence-Lewis  (1977)  present  stationary  M.A.  models 
where  the  r.v.s  X(n)  have  exp.  dist.s  •  Gaver—  Lewis  (1978)  consider  stationary 
A.R.M.A.  models  where  the  r.v.s  X(n)  have  gamma  dist.s.  Jacobs-Lewis  (1978,a,b,c, 
1983)  construct  A.R.M.A.  models  where  the  r.v.  s  X(n)  are  discrete  and  assume 
values  in  a  common  finite  set.  The  aformentioned  models  can  be  used  in  the 
various  fields  of  applied  probability  and  time  series  analysis.  In  particular*  these 
models  have  been  used  to  model  and  to  analyze  univariate  point  processes  with 
correlated  interarrival  times.  Jacobs  (1978)  uses  the  exp.  models  in  queues 
with  correlated  service  and  correlated  interarrival  times.  Other  authors  have 


concentrated  their  efforts  on  the  analysis  of  univariate  point  processes  in  the 


context  of  time  series  analysis,  for  example  Bartlett  (1963,  1966),  Brillinger 
(1972,  1978),  Jacobs  and  Lewis  (1977)  and  Lawrance  and  Lewis  (1977)  to  mention 
a  few.  More  details  concerning  the  univariate  geo.  M.A.  process  and  the  corres¬ 
ponding  point  processes  may  be  found  in  Langberg-Stoffer  (1985). 

In  this  paper  we  present  two  classes  of  finite  and  infinite  M.A.  sequences 
of  bivariate  r.ve.s.  The  first  class  has  exp.  marginals  while  the  second  class 
has  geometric  (geo. )  marginals.  Uithin  each  class  of  M.A.,  the  sequences  are 
classified  according  to  their  order  of  dependence  on  the  past.  For  the  sake  of 
clarity  we  restrict  ourselves  to  bivariate  M.A.  sequences.  However  these  models 
can  be  extended  in  a  straight  forward  way.  We  use  the  theory  of  positive  de¬ 
pendence  to  show  that  in  a  variety  of  cases  the  two  classes  of  M.A.  sequences  are 
associated  (AS.).  We  then  apply  the  association  to  obtain  bounds  for  the  bi¬ 
variate  point  processes  and  to  establish  some  moment  inequalities. 

In  Section  2  we  define  the  bivariate  exp.  and  geo.  dist.s,  which  are  the  under¬ 
lying  dist.s  of  our  two  classes,  and  present  a  variety  of  examples  of  such  dist.s. 

Further  in  Section  2  we  define  the  concept  of  association  and  present  a  variety 
of  bivariate  exp.  and  geo.  dist.s  that  are  AS.  In  Section  3  we  construct  the 
two  classes  of  M.A.S proving  that  they  have  exp,  or  geo.  marginals  and  showing 
that  if  the  underlying  dist.  is  AS.,  sols  the  related  M.A.  sequence.  Finally  in 
Section  3  we  present  the  autocovariance  matrices  for  both  classes  of  M.A.  se¬ 
quences.  In  Section  4  we  define  bivariate  point  processes  whose  interarrival 
times  are  described  by  the  bivariate  exp.  or  geo.  M.A.  processes  discussed  in 
Section  3.  We  show  in  Section  4  that  if  the  interarrival  times  are  AS.  the  corresponding 
point  processes  inherit  positive  dependence  properties.  We  then  exploit  positive 
dependence  to  obtain  bounds  on  the  joint  probabilities  of  the  point  processes. 

We  conclude  Section  4  with  moment  Inequalities  for  the  bivariate  processes  and 
their  interarrival  times. 


2 .  Preliminaries . 


In  this  section  we  present  definitions  and  prove  some  basic  results  to  be 
used  in  the  sequel.  First  we  present  a  definition  of  a  bivariate  geometric 
distribution  (B.G.D. ) . 

Definition  2.1.  Let  N(l),  N(2)  be  r.v.s  assuming  values  in  the  set  {1,2,...}. 
We  say  that  (N(l),  N(2))  has  a  B.G.D.  if  N(l),  N(2)  have  geometric  distributions. 
Examples  2.2.  (a)  Let  N  be  a  geometric  r.v.  Then  (N,N)  has  a  B.G.D. 

(b)  Let  N(l),  N(2)  be  independent  geometric  r.v.s.  Then  (N(1),N(2))  has  a  B.G.D. 

(c)  Let  N(.l)  ,N(2)  ,N(3)  be  independent  geometric  r.v.s.  Then  (min(N(l)  ,N(3) ) , 
min(N(2),N(3)))has  a  B.G.D.;  the  Esary-Marshall  (1974)  B.G.D.  (d)  Let  P(0,0), 

P(0,1),  P(1,0),  and  P(l,l)  be  in  [0,1]  s.t.  (i)  P(0,0)  +  P(0,1)  +  P(1,0)  +  P(l,l)  «  1, 
(ii)  P(0,1) +  P(l,l) ,  P(1,0) +  P(l,l)  <1,  and  let  N(1),N(2)  be  r.v.s  assuming 
values  in  the  set  {l,2,...}  determined  by: 

T  [P(l,l)]a[P(0,l)  +  P(l,l)]b"a,  b>a 

(2.3)  P(N(1)  >  a,  N(2)  >  b)  -  |  ,  a,b«l,2 . 

I  [P(l,l)  j  [P(1,0)+  P(l,l)]a_  ,  b^a 

Then  (N(1)»N(2))  has  a  B.G.D.;  the  Block  (1977)  fundamental  B.G.D.  (see  also 
Block-Paulson  C1984)).  (e)  Let  (M(J.),M(2))  have  a  B.G.D.  and  let  N(j)  » 

(N(j,l)  ,N(j  ,2)) ,  j*l,2,...,  be  an  i.i.d.  sequence  of  random  vectors  (r.ve.s), 
independent  of  (M(l)  ,M(2)),with  a  B.G.D.  Then  (^^N(j,l),  ^^N(j,2))  has  a 
B.G.D. 

In  the  following  remark  we  show  that  Examples  (Ex.s) (a) , (b)  and  (c)  in  Ex. (2.2) 

(but  not  Ce))  are  particular  cases  of  Ex.  (2.2)(d). 

Remarks  2.4.  (a)  Let  P(1,0) *  P(0,1)  ■  0.  Then  we  obtain  the  B.G.D.  introduced 

in  Ex.  2.2(a).  (b)  Let  P(l,l)»  (P(l,l)+  P(1,0))(P(1,1)+  P(0,1)).  Then  we  obtain 

the  B.G.D.  introduced  in  Ex.  2.2(b).  (c)  Let  P(l,l)>  (P(l,l) +  P(1,0)) (P(1,1)+P(0 ,1) ) 


and  let  MCI),  M(2) ,  M(3)  be  independent  geometric  r.v.3  with  parameters 
P(1,1)(P(1,1)  +  PC1.0))"1,  PC1,1)(P(1,1)  +  P(0,1))_1,  [P(l,l)]”1(P(l,D  + 
P(1,0))(P(1,1)  +  P CO, 1)), respectively.  Then  (N(1),N(2)),  given  by  (2.3),  is 
stochastically  equal  to  the  Esary-Mar shall  B.G.D.:  (min(M(l) ,M(3)) ,min(M(2) ,M(3))) . 

Next  we  present  a  definition  of  a  bivariate  exponential  (exp. )  distribution 
(B.E.D.) 

Definition  2.5.  Let  E(l),  E(2)  be  r.v.s  assuming  values  in  (O,00).  We  say 

that  (E(1),E(2))  has  a  B.E.D.  if  E(1),E(2)  have  exp.  distributions. 

Examples  2.6.  (a)  Let  E  be  an  exp.  r.v.  Then  (E,E)  has  a  B.E.D.  (b)  Let 

E(1),E(2)  be  independent  exp.  r.v.s.  Then  (E(1),E(2))  has  a  B.E.D.  (c)  Let 

E(l) ,E(2) ,E(3)  be  independent  exp.  r.v.s.  Then  (min(E(l) ,E(3)) ,  min(E(2) ,E(3) ) ) 

has  a  B.E.D.;  the  Marshal l-01kin  (1967)  B.E.D.  (d)  Let  (N(1),N(2))  have  a 

B.G.D.  and  let  E(j)-  (E(j,l)»  E(j,2)),  j»l,2,...,  be  an  i.i.d.  sequence  of 

r.ve.s,  independent  of  (N(l)  ,N(2) ),  with  a  B.E.D.  Then  ,1) ,  Ijl^E(j,2)) 

has  a  B.E.D.  (e)  Let  0  <_  a  <1.  Then  (E(1),E(2))  determined  by:  P(E(1) >  x, 

E(2)>y)  -  e-x-y-°«y  ^  x ,y  >  <*  has  a  B.E.D.;  a  Gumble  (1960)  B.E.D.  (f)  Let 

|ot |  <1.  Then  CE(1),E(2))  determined  by:  P(E(1)  <  x,  E(2)  <y)  -  (l-e~x) (l-e_y) 

(1+ae  X  y)*  x,y  >  0,  has  a  B.E.D.;  a  Gumble  (I960)  B.E.D.  (g)  Let  a_>  1.  Then 

a  a. 1/a 

(E(l) ,E(2) )  determined  by:  P(E(1)  >  x,  E(2)  >  y)  ■  e  x  “  ,  x,y  >  0,  has 

a  B.E.D.  (h)  Let  (X,Y)  be  a  r.ve.  with  continuous  marginal  distributions  F,G 
respectively.  Then  (-ln(l-F(X)) ,  -lnCl-G(Y)))  has  a  B.E.D. 

Ex.  2.6(d)  has  been  used  by  several  researchers  to  generate  bivariate  distributions 
(for  example  Arnold  (1975),  Downton  (1970),  and  Hawkes  (1972)  to  mention  few). 

In  the  following  remark  we  illustrate  how  some  of  the  B.E.D.  are  obtained  from 
Ex.  2.6 Cd) . 
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Remarks  2.7.  (a)  Let  N(l)  -  N(2)  and  let  E(j,l),  E(j,2)  be  independent 
r.v.s,  j*  1,2,...  .  Then  we  obtain  the  B.E.D.  introduced  by  Downton  (1970). 

(b)  Let  (N(l) ,  N(2))  be  as  in  Ex.  2.2(d)  and  let  E(j , 1) ,  E(j,2)  be  independent 
r.v.s,  j-  1,2,....  Then  we  obtain  the  B.E.D.  introduced  by  Hawkes  (1972)- 
Paulson  (1973).  (c)  Let  (N(1),N(2))  be  as  in  Ex.  2.2(c)  and  let  E(j,l)  - 

E  CJ » 2) ,  j*  1,2,....  Then  we  obtain  the  Marshall-Olkin  (1967)  B.E.D.  given  in 
Ex.  2.6(c)  (for  details  see  Marshall-Olkin  (1967)). 

Finally  we  present  a  concept  of  positive  dependence. 

Definition  2.8.  Let  T  »  (T(H,. . .  ,T(n).)_ ,  n*l,2,...,  be  a  multivariate  r.ve. 

We  say  that  the  r.v.s  T(l) ,. . . ,TCn)  are  associated  (AS.)  if  for  all  pairs  of 
measurable  bounded  functions  f,g;  Rn-+R  both  nondecreasing  in  each  argument 
cov(f  (T)  ,g(T) )  >0. 

Remarks  2.9.  (a)  Note  that  independent  r.v.s  are  AS. (Barlow-Proschan 

(B-P)  (1975),  Th.  2.2  p.  31).  and  that  nondecreasing  functions  of  AS.  rv.s  are 
AS .  (B-P  P3  p.  30).  Thus  the  components  of  the  r.ve.  given  in  Ex.~~2.2(c)  and  the  com¬ 
ponents  of  the  r.ve.  given  in  Ex.  2.6(c)  are  AS.  (b)  Let  (E(1),E(2))  be  as  in  Ex.  2.6(e) 
(with  a  >  0)  or  in  2.6(f)  (with  -l<a<0).  Since  P{E(1)  >  x,  E(2)>y)  <  P{E(1)  >  x} 

P{E(2)  >y}  for  x,y  >  0,  E(l),  E(2)  are  not  AS.  (d)  Let  (X,Y)  be  as  in 
Ex.  2.6(h).  If  X,Y  are  AS.  (not  AS.)  then  -In  (l-F(X)),  -ln(l-G(Y))  are 
AS.  (not  AS.)  (by  P3  of  B-P  p.  30). 

The  following  lemma  provides  sufficient  conditions  for  some  of  the  bivariate  dist.s 
presented  in  Ex.s  2.2  and  2.6  to  be  AS. 

Lonrna  2.10.  Let  Q  -  (Q(l) ,Q(2) )  be  a  r.ve.  with  components  assuming  values 
in  the  set  {1,2,...}  and  let  R(j )  -  (R(j  ,1)  ,R(j  ,2)) ,  j-1,2,...  be  an  i.i.d.  sequence 
of  r.ve.s  with  nonnegative  components  independent  of  Q.  Assume  that  Q(1),Q(2)  are 
AS.  and  that  R(l,l),  R(l,2)  are  AS.  Then  £Q(1)R(.1 ,1)  ,IQ(2)R(j  ,2)  are  AS. 
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Proof :  Let  f,g:  R  +R  be  measurable  bounded  functions  nondecreasing  in  each 
argument  and  let  U(£)  ■  J^^R(j  ,£) ,  £,  =  1,2.  First  note  that 

j-l 

cov{f(ua),U(2)),g(U(l),U(2)}  =  E{cov(f(U(l)U,(2)),g(U(l),U(2))|Q} 

+  cov{Ef (U(l) ,U(2)) |q,  EgCL(l) ,U(2)) |Q>. 

Now  Ef (U(l) ,U(2) ) Iq,  Eg(U(l) ,U(2)) )Q  are  nondecreasing  functions  of  Q(1),Q(2). 
Since  Q(1),Q(2)  are  AS. 

cov{Ef (U(l) ,U(2)) |q,  Eg(U(l),U(2))  |q)}  >_0. 

Next  let  Q  *  max  (Q(1),Q(2)),  f (U(l) ,U(2) ) |Q ,  g(U(l) ,U(2) |Q  are  nondecreasing 
functions  of  R(l,l) , . . . ,R(Q,1) ,  R(l,2) . ,R(Q,2) ,  by  P3and  Th.  2.2  of  B-P 
p.  30-31  these  r.v.s.  are  AS.  Thus: 

covtf  (U(l)  ,U(2)  |Q,  g(U(l)  ,0(2))  |Q}  >0. 

Consequently  cov{f (U(l) ,U(2)) ,  g(U(l) ,U(2) ) } ^0.  || 

Remark  2.11.  In  particular  we  conclude  from  Lemma  2.10  that:  (i)  The 
components  of  the  B.C.D.  given  in  Ex.  2.2(e)  are  AS.  provided  M(1),M(2)  are  AS. 
and  N(l,l),  N(l,2)  are  AS.  (ii)  The  components  of  the  B.E.D.  given  in  Ex.  2.6(d) 
is  AS.  provided  N(1),N(2)  are  AS.  and  E(l,l),  E(l,2)-  are  AS. 

3.  Model  Constructions. 

In  this  section  we  construct  two  classes  of  finite  and  infinite  M.A.  se¬ 
quences  of  bivariate  r.ve.s.  We  denote  the  first  class  of  sequences  by 
(X(n,m)  ■  (X(n,m,l) ,  X(n,m,2)),  n  *  0,+l,+2». . .  }  m  *  1,2, and  the  second  class 
of  sequences  by  (G(n,m)  ■  (G(n,m,l) ,G(n,m,2)) ,  n  »  0,+l,+2, . . . }  ra=  1,2,. . . We 
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show  that  each  r.ve.  X(n,m)  has  a  B.E.D.  with  a  vector  mean  that  does  not  depend 
on  n  or  m  and  that  each  r.ve.  G(n,m)  has  a  B.G.D.  with  a  vector  mean  independent 
of  n  or  m.  Within  each  class  of  sequences  the  order  of  dependence  on  the  past 
is  indicated  by  the  parameter  m.  For  each  positive  integer  m,  X(n,m)  (G(n,m)) 
depends  only  on  the  previous  m  r.ve.s:  X(n-l,m) , . . . ,X(n-m,m)  (G(n-l,m) , . . . ,G(n-m,m)) 
while  the  r.ve.  X(n,°°)  (G(n,°°) )  depends  on  all  the  preceeding  r.ve.s:  X(n-1,«>), 

X(n-2 ,«),...,  (G(n-1 ,°°) ,  G(n-2,°°) , . . . )  .  After  constructing  the  various  models  we 
present  sufficient  conditions  for  the  r.v.s.  X(n(j )  ,m,£)  (G(n(j )  ,m,£)),  £  =  1,2, 
j  =  1, . . .  ,k  to  be  AS .  (where  k  =  1,2, . . . ,  and  n(l)  <n(2)  <...<  n(k)  e  (0,+l,+2, . . .  ’). 

We  conclude  this  section  by  computing  the  autocovariance  matrices  for  the  two 
classes  of  sequences  and  present  sufficient  conditions  for  the  sequences  to  be 
stationary.  For  the  stationary  sequences  we  give  the  spectral  density  matrices. 

First  we  construct  the  "exp."  class.  Some  notation  is  needed. 

Notation  3.1.  Throughout,  n  ranges  over  the  integers  and  m,j  over  the 
positive  integers.  Let  E(tij)  =  (E(n,l) ,  E(n,2))  be  i.i.d.  r.ve.s  with  a  B.E.D. 
and  mean  vector  (A  1(1),A  *(2))  (A(1),A(2)  >0).  Let  (B(n, j ,1) ,B(n,j ,2) )  be  bi¬ 
variate  vectors  with  components  in  [0,1]  and  let  B(n,j)  be  a  2*2  diagonal  matrix 
with  8(n,j,l),  B(n,j,2)  on  the  main  diagonal.  Further  let  (I(n,j ,1) ,I(n,j ,2))  be 
independent  bivariate  r.ve.s  independent  of  all  the  r.ve.s  E(n)  s.t.  I(n,j,l), 
I(n,j,2)  are  Bernoullis  r.v.s  with  parameters  l-8(n,j,l),  1-B(n,j ,2) , respectively, 
and  let  V(n,j,q)  be  a  2*2  random  diagonal  matrix  with 

j  j 

II  I(n,£,l),  II  I(n,£,  2)  on  the  main  diagonal  qe{l,...,j}. 

£=q  £=q 

Finally ,  let  a  sum  (product)  over  an  empty  set  of  indices  be  equal  to  zero  (one) . 

We  present  now  the  class  of  "exp"  sequences.  For  m»l,2,...,  and  n  =  0,+l,+2, . . . , 


(3.2) 


X(n,m)  =  V(n,r ,l)B(n,r+l)E(n-r)  +  V(n,m,l)E(n-m) , 

and 

(3.3)  X(n,<*>)  *  ^0V(n,r,l)B(n,r+l)E(n-r). 

We  show  in  Corollary  3.8  and  Lemma  3.9  that  for  all  n,m,  X(n,m)  and  X(n,“)  have 

B.E.D.s.  Next  we  construct  the  "geometric"  class.  Some  notation  is  needed. 

Notation  3.4.  Let  p(l),p(2)  be  real  numbers  in  (0*1]  and  let  (a (n,l) ,a (n, 2) ) 

be  bivariate  vectors  s.t.  p(£)  <_  ot(n,2)  <_1,  2 ,  ■  1,2.  Further  let  N(n)  *  (N(n,l)  ,N(n,2) ) 

be  independent  r.ve.s  with  B.G.D.s  and  mean  vector  (p  ^ (l)a(n,l) ,  p  ^(2)a(n,2)) 

and  let  M(n)  =  (M(n,l) ,M(n,2))  be  i.i.d.  r.ve.s,  independent  of  all  N(n), 

-1  -2 

with  a  B.G.D.  and  mean  vector  (p  (l),p  (2)).  Finally  let  (J(n,j ,1) , J(n, j ,2) ) 

be  independent  r.ve.s  independent  of  all  previous  r.ve.s  s.t.  J(n, j ,1) , J(n, j ,2) 
are  Bernoulli  r.v.s  with  parameters  l-a(n,l),  l-a(n-2),  respectively, and  let 

j  j 

U(n,j,q)  be  a  2*2  random  diagonal  matrix  with  II  J(n,2,l),  n  J(n,2  ,2)  on  the 

2=q  2=q 

mam  diagonal  qe{l,...,j). 

We  now  present  the  class  of  "geometric"  sequences.  For  t*l,2,...,  n  *  0,+l,+2, . . . , 
let 

(3.5)  G(n,m)  =  U(n,r,l)N(n-r)  +  U(n,m+1, l)M(n-m) 

and 

(3.6)  G(n,°°  )  =  J”  U (n ,r , l)N(n-r) . 

r=0 

Next  we  show  that  X(n,m)  (G(n,m))  has  a  B.E.D.  (B.G.D. ).  The  following  lemna  is 
needed. 

Lemma  3.7.  Let  Y(n,m,q)  =  Z^=q  V(n,r+q-l,q) B(n, r+q) E(n-r-q+l)  + V(n ,m+q-l ,q ) 
E(n-m-q+l),  (W(n,m,q)  *  U(n,r+q-l,q)N(n-r-q+l) + U(n,uH-q ,q)M(n-m-q+lJ, 

n  ■  0,+l,+2 , . . . ,  ra,q  ■  1,2 .  Then  for  all  n,m,  and  q,  Y(n,m,q) (W(n,m,q ) )  has  a 


V  s 
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B.E.D.  (B.G.D.)  with  mean  vector  (X  ^(D.X  *(2))  ((p  ^(lKp  *( 2 ))). 

Proof:  We  prove  the  result  of  the  lemma  by  an  induction  argument  on  m. 

For  m«l  (m«0)  Y(n,l,q)  -  B(n,q)E(n-q+l)  +  V(n,q,q)E(n-q)  (W(n,0,q)  -  N(n-q+l)  + 
U(n,q,q)M(n-q+l) ) .  By  computing  the  characteristic  functions  of  the  components 
of  Y(n,l,q) (W(n,l,q))  one  can  verify  that  the  results  of  the  lemma  hold  for  all 
n,q.  Let  us  assume  that  the  results  of  the  lemma  hold  for  m  andalln.q.  Note 
that  Y(n,nri-l,q)  -  B(n,q)E(n-q+l)  +  V(n,q,q)[£^g  V(n,r+q,q+l)B(n,r+q+l)E(n-q-r) 

+  V(n,mfq,q+l)E(n-m-q)](W(n,nH-l,q)  -  N(n-q+l)  +U(n,q,q)  [^QU(n,r+q,q+l) 

N(n-q-r)  +  U(n,m+q+l,q+l)M(n-m-q)  ]) .  By  the  induction  assumption  the  r.ve. 
in  the  brackets  has  a  B.E.D.  (B.G.D.)  with  mean  vector  (X  *(1),X  ^(2)) 

((p  ^(l),p  "'■(Z))).  Since  this  r.ve.  is  independent  of  E(n-q+l)  (N(n-q+l))  it 
follows,  as  in  the  case  m*l  (m-0) ,  that  Y(n, art-1, q)  (W(n,nrt-l,q))  has  a  B.E.D.  (B.G.D.) 

•*>  tm 

with  mean  vector  (X  1(1),X  1(2))((p  ^(D.p  1(2)))  for  all  n,q.  || 

Note  that  X(n,m)  (G(n,m)),  given  by  (3.2)  ((3.5)),  is  equal  to  Y(n,m,l)  (W(n,m,l)). 
Thus  we  conclude  from  Lemma  3. 7  that 

Corollary  3.8.  For  all  n,m  X(n,m)  (G(n,m))  has  a  B.E.D.  (B.G.D.)  with 
mean  vector  (X_1(l) ,x"1(2)) ( (p_1(l) ,p"1(2)) ) . 

We  show  now  that  G(n,«>),  given  by  (3.6),  has  a  B.G.D.  Also  we  show  that  if  for 
all  n  and  i *  1,2 

m 

(3.9)  lim  n  (l-6(n, £.,!))  -  0 

m-*»  lml 

then  X(n,«»),  given  by  (3.3),  has  a  B.E.D. 

Lemma  3.10.  (a)  For  all  n,G(n,«)  has  a  B.G.D.  with  mean  vector  (p  ^(D.p  ^(2)). 
(b)  If  Condition  (3.9)  holds  then  for  all  n,X(n,®)  has  a  B.E.D.  with  mean  vector 
(X'1(1),X"1(2)). 


Proof:  Let  n  be  a  positive  Integer.  Since  lim  (l-a(n,i))m  <_  lim(l-p(i) )m 


particular  X(n,m)  — 2->  x(n,»)  and  G(n,m)  G(n,®).  Consequently  the  results 

m  nr+°°  _  ,  nr*®  - 

of  the  lemma  follow  from  Corollary  3.8.  | | 

Note  that  for  (3.9)  to  hold  it  suffices  that  for  all  n  and  i  *1,2, 
inf (B(n,4,i) ,  4  ■  1,2  ,...}>  0. 

Next  we  investigate  some  dependency  aspects  of  both  classes. 

Remarks  3.11.  (a)  For  a  fixed  m,  the  sequences  (X(n,m),  n  -  0,+l,+2,. . . }, 
{G(n,m),  n * 0,+l,+2, . . . }  are  m-dependent  (i.e.  for  n(l),  n(2)  integers  s.t. 
|n(l)-n(2)|  >m  the  r.ve.s  X(n(l),m),  X(n(2),m)  (G(n(l),m),  G(n(2),m))  are  inde¬ 
pendent).  (b)  Clearly  if  we  choose  m  to  be  a  function,  say  of  n  (v(n)  e  {1,2,.. 
for  all  n)  then  the  dependency  of  X(n,  ^ (n) )  (G(n,iKn)))  on  the  past  varies  with 
n.  (c)  It  is  easy  to  see  that  for  all  n  the  r.ve.  X(n,«)  (G(n,“))  depends  on  all 
the  preceeding  r.ve.s  X(q,«),  -  “  <  q  <  n  (G(q ,<*>),  -  ®<q<n). 

We  now  investigate  one  positive  dependence  aspect  of  both  classes. 

Lemma  3.12.  Suppose  that  E(l,l) ,E(1,2)  are  AS.’  Then  for  all  positive  In¬ 
tegers  m,k,  and  all  integers  n(l)  <n(2)  <...  <n(k)  the  r.v.s  X(n(j),m,4), 

4*1,2,  j  *  1, . . . ,k  are  AS. 

Proof :  By  Th.  2.2  p.  31  and  p.  30  of  R-P,  the  r.v.s.  E(r,4),  I(n(j),q,4), 
4*1,2,  q*l,...,m,  j-l,...,k,  r -n(l)j...  ,n(k)  are  AS.  Since  the  r.v.s 
X(n(j),m,4),  4-1,2,  j-l,...,k  are  nondecreasing  functions  of  the  previous  collec¬ 
tion  of  AS.  r.v.s  the  result  of  the  lemma  follows  by  P^  p.  30  of  B-P.  | | 

In  a  similar  way  one  can  prove  the  following  lemma. 

Lemma  3.13.  Suppose  that  M(l,l),  M(l,2)  are  AS.  and  that  for  all  n, 

N(n,l),  N(n,2)  are  AS.  Then  for  all  positive  integers  m,k,  and  all  integers 
n(l)  <  n(2)  <  . . .  <  n(k)  the  r.v.s  G(n(j)  ,m,4)  ,  4-1,2,  j-l,...,k  are  AS. 

Next  we  prove  similar  results  for  the  sequences  {X(n,«),  n  -  0,+l,+2,. . . }  and 
(G(n,®),  n  -  0,+l,+2,. . . }. 

Lemma  3.14.  (a)  Let  us  assume  that  M(l,l),  M(l,2)  are  AS.  and  that  for  all 
n,  N(n,l),  N(n,2)  are  AS.  Then  for  all  positive  integers  k  and  all  integers 


n(l)  <n(2)  <  ...  <n(k),  the  r.v.s  G(n(j),®,£),  £-  1,2,  j  ■  l,...,k  ere  AS.  (b)  Let 
us  assume  that  E(l,l),  E(l,2)  are  AS  and  that  Condition  (3.9)  holds.  Then  for 
all  positive  integers  k  and  all  integers  n(l)  <n(2)  <  ...  <n(k),  the  r.v.s 
X(n(j),»,£),  £-1,2,  j-l,...,k  are  AS. 

Proof :  By  similar  arguments  to  the  ones  given  in  the  proof  of  Lemma  3.10 
we  conclude  that  the  two  sequences  of  r.ve.s:  (G(n(l) ,a,l) ,  G(n(l) ,m,2) , . . . , 

G(n(k),m,l),  G(n(k) ,m,2))  and  (X(n(l) ,m,l) ,  X(n(l),m,2) . X(n(k),a,l), 

X(n(k)  ,m,2)) ,  m-1, 2,...,  converge  in  distribution  as  m-*®  to  the  r.ve.s: 

(G(n(l) ,®,  1),  G(n(l) ,®,2) , . . . ,G(n(k) ,®,1) ,  G(n(k),-,2)  and  (X(n(l) ,«,1) , 
X(n(l),®,2),...,X(n(k),®,l),  X(n(k) ,®, 2)),  respectively.  By  Lemma  3.12  the  r.v.s 
X(n(j),m,£),  £-1,2,  j-l,...,k,  are  AS.  for  all  m  and  by  Lenma  3.13  the  r.v.s 
G(n(j)  ,m,£)  ,£  -  1,2,  j-l,...,k,  are  AS.  for  all  m.  Consequently  the  results  of 
the  lemma  follow  by  of  Esary-Proschan-Ualkup  (1967).  || 

Next  we  confute  the  autocovariance  matrices  for  both  classes  of  sequences. 
Some  notation  is  needed. 

Notation  3.15.  Let  r(e),  T(g),  and  T(g,n)  be  the  covariance  matrices  of 

E(1),'M'(1)  and  N(n)  respectively.  For  £-0,1,...,  let  r(e,n,£,m)  and  r(g,n,£,m) 

be  the  auto covariance  matrices  of  X(n,m),  and  G(n,m),  respectively,  n,£®0,+l,+2 

«•  "» 

..., m-1, 2, ...,«.  Further  let  A(n,;j)  be  a  2x2  diagonal  matrix  with 
(l-a(n,l))^,  (l-a(n,2))^  on  the  diagonal,  let  I  be  the  2x2  identity  matrix,  and 
X  the  indicator  function.  By  some  simple  calculations  we  obtain  for  n-0,+1,+2, 
m-1, 2,...,®,  and  £-1,2,...  (but  not  zero),  that 

(3.16)  r(e,n,£,m)  - 

r  r+£ 

r^0  o(n,r+l)  [  II(I-B(n,j) )  ]r  (e)  [  n  (I-B(n+£,j)]B(n+£,r+£-l) 
j-1  J-l 

m-£  m 

+B(n,m-£+l)[  H  Cl-B(n,j )) ]r(e) [  n  (I-B(n+£,j)) ] . 
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We  may  obtain  the  off-diagonal  elements  of  r(e,n,0,m)  from  Equation  (3.16)  by 

setting  B(n,nH-l)  ■  I  and  £  •  0.  The  diagonal  elements  are  the  variances  of 

-2  -2 

X(n,m,l)  and  X(n,m, 2),  namely,!  (1)  and  X  (2),  respectively. 

In  a  similar  way  we  obtain  for  n  ■  0,+l,+2, . . . ,  m  ■  1,2, . . .  ,<»,  and  £«1,2,...  (but 
not  zero) ,  that 

(3.17)  r(gtnf£fm)  »  lrmQ  A(nf r)T (g,n-r)A(n+fc,r+£) 

+  A(n,mfl)r(g)A(n,mfl). 

We  may  obtain  the  off  diagonal  elements  of  r(g,n,0,m)  from  Equation  (3.17)  by 
setting  £-0.  The  diagonal  elements  are  the  variances  of  G(n,m,l),  G(n,m, 2),  namely, 
(l-p(l))p-2(l)  and  (l-p(2))p“2(2) ,  respectively. 

Finally  we  give  sufficient  conditions  for  the  exponential  and  geometric  se- 
quences  to  be  stationary  and  present  their  spectral  density  matrices. 

First  we  address  ourselves  to  the  exponential  case.  Let  us  assume  that 
B(n, j ,1) ,  B(n,  j ,2) ,  given  in  Notation  3.1,  do  not  depend  on  n.  Then  clearly  the 
bivariate  exponential  sequences  given  by  (3.2)  and  (3.3)  are  stationary.  Let 
us  denote  the  autocovariance  function  of  the  stationary  X(n,m)  by  r(e,£,m). 

m 

We  obtain  T(e,£,m)  from  Equation  (3.16)  by  simply  suppressing  the  index  correspond¬ 
ing  to  n  in  Equation  (.3. 16).  Note  that  for  m<®  (e»A,®)||  1®  a  finite  sum 

and  is  thus  finite,  and  for  m«<*>,  J  ||r  (e,£,m)||  <»by  Cons  it  ion  (3.9).  Consequently 

the  2x2  spectral  density  matrix,  f(e,oi,m),  of  the  stationary  process  X  fa, m)  is  given 
by: 

(3.18)  f(e,w,m)  -  -jjj-  <*>€  I“1T*1T3* 

m  ■  1, . . . ,®, 


with  inverse  relationship 


(3.19) 


r(e,£,m)  m  I  f(e,u,m)e“  ”vdo>,  £  *  0,+l,+2, . . . , 

J  —  TT 


d  *  1,2, ... ,c 


(where  r(e,-£,m)  is  Che  transpose  of  r(e,£,m)). 

A  sufficient  condition  for  the  geometric  sequences,  given  by  Equations  (3.5)  and 
(3.6),  to  be  stationary  is  that  a(n,l),  a(n,2),  given  in  Notation  3.4  are 
equal  to  a(l),  a(2)  respectively  for  all  n.  Let  us  denote  the  autocovariance 
function  of  the  stationary  sequence  G(n,m)  by  r(g,£,m).  We  obtain  r(g,£,m)  from 
Equation  (3.17)  by  suppressing  the  index  corresponding  to  n  in  Equation  (.3.17).  Si 
a(l)  >_p(l)  and  a(2)^p(2),  mll r  (g,£,m)||  <  «  for  m  «  1,2, . . .  Consequently  the 
2x2  spectral  density  matrix,  f(g,u>,m),  of  the  stationary  process  G(n,m)  is  given  by 


(3.20) 


f(g,o),m)  -  Z1»_a,I'(g,i,m)e  *  ,  u>  e  [-it,ir]. 


m  -  1,2,...,», 


with  inverse  relationship 


<3*21)  r(g,£,m)  -j  f(g,w,m)e,'“*’du  ,  £«0,+l,+2,. 

m  -  1,2,...,», 

(where  r(g,-£,m)  is  the  transpose  of  r(g,£,m)). 


4.  Inequalities. 

For  m-l,2,...,»,  let  {R(s,t,m) «  (R(s,m,l),  R(t,m,2)),  s,t  >0}, 

R(0,0,m)  *  (0,0) ,  be  a  bivariate  point  process  with  interarrival  times  equal  to 
X(l,m),  X(2,m) ,. . . ,  given  by  Equation  (3.2)  or  (3.3).  Furthermore  let 
{S(a,b,m)  •  (S(a,m,l),  S(b,m,2)),  a,b  ■  0,1,. . . },  S(0,0,m)  -  (0,0),  be  a  bivariate 
point  process  with  interarrival  times  equal  to  G(l,m) ,G(2,m) , . . . ,  given  by 
Equation  (3.5)  or  (3.6).  He  show  that  if  the  interarrival  times  of  . the  process 
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Now  we  show  that  if  the  interarrival  tines  of  R(S)  are  AS.,  then  R(S) 
inherits  positive  dependency  properties. 

Lemma  4.6.  Let  m  -  1,2, . . . ,«,  and  let  us  assume  that  for  q-1,2,...,  the 
r.v.s  {X(j,m,£),  j  - 1,. . .q,  1*1,2}  are  AS.  Then  for  all  positive  real  numbers 
s » •  •  •  »s^  and  -  1,2  the  r.ve.s  (R(s^,m,ij),  j-l,...,q)  are  P.U.O.D. 

and  P.L.O.D. 

Proof :  Let  n^,...,n^  be  positive  Integers,  let  f  * x  n 

(Zr^1X(r,m,i^)<  Sj) 

and  let  g  -  x  _  ,  J-l,...,q.  The  functions  f.  ,...,f  (g,  ,...,g  ) 

J  n.  1  q  1  q 

(  Ir-lX(r,m,V“  8^) 

are  nonincreasing  (nondecreasing)  functions  of  AS.  r.v.s.  By  B-P  P^  p.  30 

flt..*,f  (g.,...,g  )  are  AS.  r.v.s.  Consequently  P{R(s  ,m,i  )  >  n  ,  j»l,...,q} 
cj  *  q  ^  J  .Jj 

-E  n  f  >  nEf  -  II  P{R(s  ,m,i  )  >n.}  and  P{R(s,,m,£.)£n.,  j -l,...,q}  » 
j-1  J  s«l  J  j-1  J  J  J  3  3  3 

q  q  q 

EH  g  >  HEg  -  II  P{R(s  ,m,£  )  <n.}.  II 
j-1  3  j-1  J  j-1  3  3  3 

Lemma  4.7.  Let  m-  1,2, ... and  let  us  assume  that  for  q«lf2f...,  the 

r.v.s  {G(j,mf£)#  j  al,..«,ql  t»lf2}  are  AS.  Then  for  all  positive  integers 

t  •  •  •  •  ^  m  the  r.v.s  (ni  >  j  *  lj  • .  •  j(j}  are  AS . 

h  h  n  J  J 

Proof :  Note  that  S(nj,m,lj)  -  J^jGCr.m,^),  j-l,...,q.  Thus  the 
SCnj.m.i^'s  are  nondecreasing  functions  of  AS.  r.v.s  and  hence  are  AS.  || 

Using  these  two  last  results  one  can  obtain  a  variety  of  probability  inequalities. 
For  more  details  see  Tong  (1980). 

In  Lemmas  4.6  and  4.7  we  bounded  from  below  joint  probabilities  of  R(S)  by 
a  product  of  marginal  probabilities  of  R(»,m,£)  (S(»,m,£)).  For  the  stationary 
models  we  bound  from  below  some  of  these  marginal  probabilities  by  Poisson  and 
negative  binomial  probabilities. 


v\  .y yyyyy  y  yy  •  yyyy.y 
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Lema  4.8.  Let  us  assume  that  S(n,l,l),  0(n,l,2)  are  equal, respectively,  to 

0(1)  ,0(2)  for  all  n.  Let  m  -  1,2,. . . £-1,2,  s>0  and  let  L(£,s)  be  a  Poisson 

_1  8 

r.v.  with  oiean  A(£)0  (£)s.  Then  R(s,m,£)  L(£,s). 

Proof:  From  Equation  (3.2)  or  (3.3)  we  obtain  that  X(q,m,£) ^0(£)E(q,£) , 

q«l,2 .  Now  for  r  *  1,2, . . .  P{R(s,m,£)  < r}  -  P{  £q-1X(q,m,£)  >  s}  1 

P{Yr  18(£)E(q,£)  >s}  -  P{L(£)  <r}.  || 

q*l 

Lemma  4.9.  Let  us  assume  that  ct(n,l),  a(n,2)  are  equal  to  a(l),  a(2), 

respectively  for  all  n.  Let  m-l,2,...,°°,  £*1,2,  r*l,2,...,  and  let 

Q(£,r)  be  a  negative  binomial  r.v.  with  parameters  (r,p(£)a  *(£)).  Then  S(r,m,£) 
s 

1  Q(£,r) . 

Proof :  From  Equation  (3.5)  or  (3.6)  we  obtain  that  G(q,m,£)  ^N(q,£) , 

q-1,2 .  Now  for  a- 1,2,...  P{S(r,m,£)  >  r+a}  -  P{[^1G(q,m,£)  >_r+a}  >_ 

P{J-r  N(q,£)  >  r+a)  >  P{Q(£,r)  >  r+a).  || 

q*X 

Using  the  AS.  of  E(l)  (G(l))  we  obtain  the  following  "residual"  inequalities. 

«* 

Lemma  4.10.  Let  the  assumption  and  notation  of  Lemma  4.8  hold  and  assume 
that  E(l,l),  E(l,2)  are  AS.  Then  for  ^  >  s^  t2-s2*  ni'n2’  re 
nx>_r,  n2>.r 

P{R(t£,m,£)  <.n£,  £ « 1,2 |R(s£,m,£)  <r,£-l,2} 

2 

>_  n  P{L(£,t£»  9£)  <_  n^-r). 

£*1 

Proof:  Since  X(q,m,£)  £0(£)E(q,£) ,  q*l,2,...,  £*1,2,  we  obtain  that 

n£+l 

P{R(t£,m,£)  <.n£,  R(s£,m,£)<r,  £  - 1,2}  -  P(][  j.  X(q,m,£)  >  t£,  ^q-1X(q,m,£)  >  s£, 

n  +1  r  n  +1 

£-1,2}  >P{fq;r+1X(q,«,£)>tr8t,  Jq-iX(q*B,£)  >«£»  **^2}  -  P{£q-r+l6(il) 

E(q,£)  >  t£-s£,  £*-jX(q,m,£)  >  s£,  £-1,2).  Note  that  Iqfr+1S(£)E(q,£)  ,  £-1,2 

are  nondecreasing  functions  of  the  AS.  rv.s  (E(q,£),  q-r+1,...,  Hj+^+I,  £*1,2} 
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and  thus  are  AS.,  and  that  they  are  Independent  of  £r^X(q,m,  £.)-,£  -  1»2.  Consequently 

q*  2 

PtR^.m.i)  <_  n^,  R(sltm,l)  <  r,  l  «  1,2}  >_  P(£*  ^(q.m,*)  >  s  ,  £«  1,2}  n 


n.+l  r  2 

P{^q-t+lSa)E(q»Jl)  >  V8l}  ”  P{VlX(q»m»)°  >  S£*  t-1,2}  n  P{L(t,tt-st)  <nt-r}. 

Now  the  result  of  the  lemma  follows.  | | 

In  a  similar  way  one  can  show  the  following  result. 

Lemma  4.11.  Let  the  assumption  and  notation  of  Lemma  4.9  hold  and  assume 

that  N(l,l) ,  N(l,2)  are  AS.  Then  for  n^,  n2>  k,  ,  r2>  a^,  a2  positive  integers 

s.t.  ri»r2>k»  n  >a^,  n2  >  a2  both  conditional  probabilities  P{S(r^,m,J.)  > 

n«.  +  ri*  lm  >  aJi+k*  a“l»2},  and  PfSCr^.m.i)  >  n£ +  r£,  l  -  l,2|S(k,m,i)  - 

a^+k,  £-1,2}  are  not  smaller  than  n  P{Q(Jl,rJl-k)^  n^  -  a^  +  r^  -  k}. 

*'  £-1 

Finally  we  address  ourselves  to  some  moment  inequalities. 

Lemma  4.12.  Let  the  assumptions  and  notation  of  Lemma  4.6  hold  and  let 

q  k  q  k 

k  , . . .  ,k  be  positive  integers.  Then  E  II  [R(s  ,m,£  )]  3  >_  II  E[R(s  ,m,£  )]  J. 
i  q  J  J  j.!  J  J 

Proof :  The  result  of  the  lemma  follows  from  Lemma  4.6  and  Inequality  (4.3) 


(with  f^(x)  ■  x  x  >_0,  j  -  l,...,q).  |  | 

Lemma  4.13.  Let  us  assume  the  assumptions  and  notation  of  Lemma  4.8  hold 

Ic  k 

and  let  k  be  a  positive  integer.  Then  ER(s,m,£)  <_  EL(£,s)  . 

Proof:  The  result  of  the  lemna  follows  by  Lemma  4.8  and  Remark  4.5  (with 
f  (x)  »xk,  x  >  OL 

In  a  similar  way  we  can  prove  the  following  result. 

Lemma  4.14.  Let  us  assume  that  the  assumptions  and  notation  of  Lemmas  4.7 

q  k1 

and  4.9  hold  and  let  k,,...,k  be  positive  integers.  Then  E  II  [S(n.,m,£,)]  J 

1  q  j-1  J  j 

q  k1 
>  n  E[Q(£  ,n  )]  J. 

j-1  J  J 

Lemma  4. 15.  Let  us  assume  that  the  assumption  and  notation  of  Lemma  4.6 

9  k4 

hold  and  let  k,,...,k  be  positive  integers.  Then  E  n  [X(j,m,£.)l  J  > 

1  q  j-1  j 


[\ 
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q  -K 

n  {k.![x(i.)]  3}. 

j-i  3  J 

Proof ;  The  result  follows  by  Corollary  3.8  and  the  association  of  the  r.v.s 
j  ■  1  * •  •  •  »q ) •  || 

Similarly  we  obtain  the  following. 


hold 


Lemma  4.16.  Let  us  assume  that  the  assumption  and  notation  of  Lemma  4.7 

q  k 

and  let  k..,...,k  be  positive  integers.  Then  E  II  [G(j  ,m,£.)  ]  J  >_ 
q  1-1  J 


q  k 

n  E[G(l,ra,£  ) ]  K 
j-1  J 


-1, 


Note  that  G(l,m,£^)  is  by  Corollary  3.8*  a  geometric  r.v.  with  mean  p  (£^), 

j  -  i»***  »q» 
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